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Abstract —Machine type communications (MTC) in 3G/4G 
networks is getting more attention recently due to bursty nature 
of traffic characteristics in contrast to Poisson type H2H traffic. A 
large number of methods have been suggested in the literature. In 
this paper we give a mathematical model on performance analysis 
of Disjoint Allocation (DA) and Joint Allocation (JA) methods for 
allocating preambles to M2M and H2H users. In an earlier work 
we had investigated the performance of two possible splitting 
preamble methods on collision probability and energy reduction 
for MTC subscribers. In this paper we develop a stochastic 
model for JA/DA method in RACH procedure using a Kth order 
Markov chain approach and carry out performance analysis in 
terms of collision probability, access success probability, average 
access delay, statistics of preamble transmissions and statistics 
of access delay for M2M users. The optimal number of reserved 
preamble set is derived based on the given success access delay 
threshold. Numerical results verifies the simulation results and 
demonstrates that this model can estimate the performance of 
the JA/DA method accurately. 

I. Introduction 

Long term evolution (LTE)/LTE-Advanced and UMTS are 
the radio access network (RAN) protocols for 4G and 3G 
standards respectively. The study of RAN improvements for 
Machine Type-Communications is still under investigation by 
the 3rd Generation Partnership Project (3GPP). This is due to 
the fact that machine service requirements are quite different 
compared to H2H services. For instance, in smart-grid net¬ 
works, Advanced Metering Infrastructure (AMI) requires low 
throughput and can tolerate high latency. Whereas, Automated 
Demand Response (ADR) requires medium throughput and 
medium to low latency and Feeder Automation (FA) requires 
high throughput and low latency. Providing the promised 
QoS for MTC subscribers under a massive deployment of 
the number of machine devices in cellular areas is what is 
known in 3GPP as RAN-overload control method. By giving 
priority to the 4G cellular networks, the service providers 
across the globe are looking for solutions which promises to 
support a wide range of applications. Towards that end, it is 
important to propose and study protocols by simulation as well 
as mathematical analysis. A variety of RAN-overload control 
methods is suggested by 3GPP m such as, extended access 
barring (EAB), separate PRACH (physical random access 
channel) resources for MTC, pull based schemes, MTC group 
paging and MTC specihc backoff scheme. No related research 
shows the effect of the separate allocating RACH resources on 
the LTE network. In this paper, the effect of preamble splitting 


between H2H and M2M users is studied with emphasis on 
modeling and performance analysis. 

Congestion or overload in RAN is in associated with a huge 
number of MTC devices attempting to establish a connection 
to the network radio component (eNB) specially in PRACH Q. 
To manage overload scenarios under the bursty arrival of the 
machine devices, 3GPP suggested splitting of the preamble 
set between H2H and M2M subscribers HI. We call it the 
Disjoint Allocation (DA) method. The splitting parameter will 
decide the performance of H2H and M2M users. However, the 
M2M subscribers will face more degradation in performance 
as a smaller number of preambles are allocated to them. 
Since the frequency of the H2H call arrival is generally much 
less compared to M2M arrival rates with bursty nature of 
arrivals we proposed and studied an enhanced version of the 
DA method where the H2H subscribers could share their 
allocated preambles with the MTC devices El. We call this 
approach. Joint Allocation (JA) method. In |[3|, we have shown 
through analysis and simulations that the JA method exhibits 
better performance in terms of the collision reduction and 
results in good energy saving for MTC devices. However, we 
mainly focused on the energy reduction for machine devices 
and not on their latency constrains. In this paper we have 
given a mathematical analysis of the RACH process under 
JA/DA method and an optimum preamble allocation to H2H 
subscribers subject to providing the required access delay for 
M2M applications. 

In order to get a better understanding of simulations and ap¬ 
proximations results of the RAN-overload control algorithms, 
recently there have been interesting analytical models in liter¬ 
ature to investigate a couple of them. We can categorize the 
RAN-overload control problems into two groups; a) resource 
allocation on PRACH and PUSCH as they share fixed number 
of RBs, where wasting any of the RBs imposes unused data 
channel on PUSCH b) only data channel allocation on PUSCH. 
In this paper our focus is on case (a), in which, we assume 
that there are proposed MAC procedures like ||9l that allows the 
eNB to be informed about the collided preambles in very hrst 
steps to prevent it allocating PUSCH RBs to them. Papers a, 
0 use an iterative model to evaluate the performance of the 
group paging algorithm in 3GPP with bursty arrival of MTC 
devices. Using mathematical models 0, 0 hnd the optimum 


'This is known as LTE upiink connection process, in which devices in 
RRC-IDLE status contend to win a data channei from eNB and turn to RRC- 
CONNECTED status O. 



group sizes and allocated radio resources for the desired QoS 
requirements of MTC devices. In ||6l author proves that steady 
state of the RACH process depends on the number of the 
permitted re-transmission times under different loads of the 
RACH arrival. Q uses a Markov chain model to effectively 
analyze the interactions in RACH channel. Analysis of the 
throughput under Poisson load in RACH procedure to examine 
the impact of different RA configurations can be found in 
ifSl . However, mathematical analysis, when the arrival of H2H 
and M2M calls following Typel and Type2 distributions as 
defined in 3GPP [I], or an analytical model to investigate the 
performance of PRACH resource allocation algorithms appears 
to be not available. In this letter, we provide a stochastic model 
and analysis for RACH procedure based on the RAN-overload 
methods DA & JA introduced in |[3l, which will be explained 
in section HVl 

The remainder of this paper is organized as follows. In section 
El we explain a preparation on the random access procedure in 
the LTE system and the related works. Section Hill elaborates 
the system model and the proposed stochastic model. Section 
HVl gives a brief introduction on the proposed JA and DA 
algorithms and analyzes the model to evaluate the performance 
metrics parameters, and we derive conditions for the optimal 
allocated number of preambles for providing the required 
access success time for machine type applications. In Section 
[Vl we present numerical results and the paper is concluded in 
Section |VT] 

H. Random Access Procedure and Related Works 
A. Random Access Procedure in LTE System 

In this subsection we provide a brief description on RACH 
procedure (RA), which proceed on PRACH 12 ■ The RA ar¬ 
chitecture of many wireless standards including IEEE 802.11*, 
IEEE 802.16 (WIMAX), UMTS and LTE is similar to a slotted 
ALOHA system. The only main difference between LTE and 
other convectional standards, is that here we are dealing with 
a multi-channel slotted ALOHA. Each RA-TS (time slot) in 
PRACH consists of multiple preambles, which are a set of non¬ 
overlapping slots (see Pigure[T]|- The RA procedure in LTE is 
completed by the following four steps. 

• Preamble transmission (stepl); A UE (user equip¬ 
ment/subscriber) randomly selects a preamble out of 
all available preambles with equal probability and 
transmits this to the eNB. 

• Random access response (step2: RAR-Msg); Upon 
detecting preambles, the eNB sends the RAR-Msg to 
the UE, which is the crucial step for granting an uplink 
resource in next step. 

• Connection setup request message (step 3); When 
a UE receives its RAR-Msg, by using the initially 
granted uplink resource from step2 the UE can send 
the connection setup request message. 

• Connection setup response message (step 4); If the 
eNB successfully receives the connection setup re¬ 
quest message in step 3, the eNB sends the acknowl¬ 
edgement message to the UE. Upon receiving this 
ACK message The RA procedure is completed. 


In Eigure[Tl we present an example of RA procedure described 
above. In this example, eight UEs arrive at {i — l)th RA-TS 
to send preamble to the eNB individually. Among them those 
preambles chosen by multiple UEs in this RA-TS (preamble #2 
& #4) are considered to be collided, and only chosen preamble 
#1 and #7 is assumed to be succesfl Those who send preamble 
without collision receive the RAR-Msg (step2) and the back- 
logged UEs (collided) should reattempt the procedure (we will 
describe the strategy of the reattempting RA procedure with 
details in section imii. 

B. Related Works 

Since the status of the RA-TSs in RA procedure can 
randomly change according to the independent arrival of 
contending subscribers, we analyze the stochastic behaviour 
of the RA procedure by design a Markov model. According 
to the literature survey we can classify the Markov analytic 
works base on the state space as follows. 

• States can represent the number of backlogged (col¬ 
lided) users states = {0,1,2, ...Njje}, in which Njje 
is the number of all users playing a contention role in 
RA process. The states can step forward or backward 
as the number of backlogged UEs can face collision 
or success respectfully lfT3ll . lfT4ll . 

• States represent the number of preamble transmission 
states = {0,1,..lU}. In which W is the maximum 
number of allowed preamble transmission a UE can 
try before it gets connected to the eNB. Otherwise, UE 
is block to participate in the contention processPH. 

• States in this approach represent the number of possi¬ 
ble RA-TSs during the experiment a UE can encounter 
states = {1, 2, .A, ..rj, S, d}. In which p is the maxi¬ 
mum RA-TS in the experiment, S is the success state 
that stores all successful UEs received their success 
RAR Msg., and state d includes the dropped UEs who 
exceeded their W number of preamble transmission 
and are withdrawn from the process ifTSl . Here, in 
the referenced article, author has either simplified the 
states’ population model by omitting the impact of 
new arrivals and number of retransmission, or the 
probability of transmission is not being considered as 
a function of channel resources. 

In this paper we introduce a order Markov chain (K-Mc) 
an enhanced version of model in na to estimate the average 
access success delay for M2M with separate preamble RAN- 
overload control method. In next section we will discuss our 
proposed K-Mc model in details. 

HI. System Model 

We consider a single eNB at the center of a cellular area 
in LTE network, in which Nh2h and Nmtc number of H2H 
and M2M subscribers are contending together. Before we start 
the analysis model, it is worth to discuss which mathematical 

^We note that there are physical layer factors such as different levels of 
transmission power among the UEs that can influence the probability of the 
success/collision of RA request. However, we focus on the behaviour of MAC 
layer, and thus the properties of the PHY layer is beyond the scope of this 
study. 
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model can fit on the MTC calls? 

In typical M2M calls, the uplink traffic dominates the downlink 
in comparison to the H2H calls (for instance sending regu¬ 
lar/irregular reports). Machines pars more number of messages 
with short length which eventually this trend generates a bulk 
amount of attempts in a cellular areas. That is the reason the 
standard Poisson model does not meet all the characteristics 
of the M2M calls. In this letter, following 3GPP H] we use 
the traffic type2, a Beta pdf for MTCs. Figure [1] illustrates 
the time evolution of the RA procedure in uplink PRACH. 
Time is divided in 1ms slots representing the full duration of 
a sub-frame (s.f) in LTE-A. RA-TS are labeled as 1 < i < 77 
. We assume that the total number of Nh 2 h + Nmtc UEs 
activate during ry time. When a UE fails to transmit its 
< 71 < W) RA request (unsuccessful preamble) at time 
slot {i — l)th it repeats the RA procedure and reattempts in 
future upcoming RA-TS {i,i + l,i + 2...). Here n = 1 refers 
a new arrival UE who is attempting the RA procedure as first 
time and n = W refers to a failed UE for W —1 times, which 
is reattempting the RA procedure for the last permitted time. 
Therefore arrival in each RA-TS can be a composition of new 
and backlogged RA arrival of H2H and M2M calls. 

If there is a new arrival (n = 1) UE departing between two 
RA-TS it waits and attempts its preamble transmission on 
immediate next RA-TS. Eor instance, in this example seven 
new arrival UEs who departure between — and ith RA- 
TSs will attempt their hrst preamble transmission on ith RA¬ 
TS. After sending the preamble message UE should wait upto 
maximum TjiAR + WiiAR time (RAR delay response and RAR 
window) to receive the RAR-Msg. In this example two UEs 
who contending in RA-TS and picked the preamble #7 

& #1 will receive their RAR-Msg after 7ms. Upon receiving 
the RAR-Msg UE is considered to be success and its total 
delay is the variance + Trar + Wrar) - 4- In which, 
tk, 4th resemble the time of UE’s initial arrival in kth TS 
(note that it can be a RA-TS or another uplink TS), and the 
successful jth RA-TS attempt. The rest of UEs who picked 
up a same preamble won’t be able to receive their RAR-Msg 
are failed. Failed UEs after Trar + Wrar time will start 
back-off procedure, by choosing uniformly a random number 
between {0-Wbo\- UE starts counting down, upon its timer is 
expired, it starts reattempting nth RA procedure (2 < n < W) 
in upcoming RA-TS. Details of all notations used in this 
paper is given in Table |I] This process continues till We 
are interested in estimating the average access success delay 
for UEs which is the interval between the hrst RA attempt 


and the RA procedure completion that will be explained in 
following section^ A K-Mc model is introduced in following 
subsection. 

A. Framework/Problem formulation 

We are interested in hnding a probability for collision 
and success transmission for each RA-TS that leads us in 
computing the expected access success delay and other per¬ 
formance metrics through analysis. By considering this fact 
that the number of contending users (M2M or H2H) per 
RA-TS is stochastic by nature, we denote them with the 
random variable, ki) as depicted in Eigure|2] In which, 

zl"\t,ki),2 < n < W is the expected number of arrival 
at given time t in RA-TS after (n — 1)*^ preamble re¬ 
transmission in past t — ki time, and its transactions can be 
modeled by a order Markov chain which we are going to 
discuss in this section. 

As it is deduced from the Figme^zj^^^ (t, ki) is determined by 
the composition of the number of new arrivals Z^^'^ (t) {n = 1) 
and number of failed users z\'^\t,ki), 2 < n < W who 
come from the i — l‘^,f — 2*^...* — previous RA-TSs. 
The cross symbols in this diagram represent the state spac43 
S = {Zi{t,ki),Z2{t,k2)-.Zi{t,ki),...Zj^{t,kr,),Zs{t)}, 

where 1 < i < rj, Zsit) is the success state that consist 
of successful UEs at time t and the sequences of time 
ti < t 2 < < T are considered to be 1ms and T is the 

latest duration of the users activation during the experiment. 
Zi{t,ki) is the expected number of users arrival at time t in 
jjj total, which is dehned as 
w 

Z,{t,h) = Y,z\^\t,h). ( 1 ) 

n—1 

To clarify the k step transmission in our model let us recap 
the back-off technique in RA procedure. Let’s assume at 
time t, number of new arrivals departure in state 

1. without contributing the failed UEs from previous states, 
these population of users will contend together on M = 

^One should note that in this paper the delay of PHY layer such as 
processing time from eNB or UE end and contention resolution time-window 
is omitted. Also as discussed earlier we have assumed that eNB can recognize 
about the fail/success status of preamble requests at Msg2, therefore the the 
delay of receiving Msg3-Msg4 is not considered in our computation. 

^Sate space is resembling the possible RA-TSs in RA procedure, therefore 
we use them interchangeably. 
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Fig. 2. A K-Mc demonstration for RACH stochastic process. 


{mi,7712, ■■■M} number of preambles which are available in 
this state. The second order failed UEs who are failed with 
their first preamble transmission will re-attend their luck in 
upcoming states (1 + 1,1 + 2 ...), and if their transmission 
failed again, they will reattempt their third luck and so on. 
UEs may continue this process of re-transmission up to W 
times. The conditional probability that the stochastic RACH 
state (RA-TS) j at time t — m will switch to state i at time t 
is 


Pjj = {X{t) = E[Zi{t,ki)]\X{t-m) =E[Zi_jn{t,ki-m,)], 
m = 1,2, ..kf, 2 < i < 7]} > 0, 

ki 

E = 1- (2) 

j=i-l 

The population of arrival in each future states depends on 
5s.f the period between two sequential RA-TS and the size 
of the back-off window Wbo- For instance, the maximum 
number of past possible steps from state that can associate 
in backlogged aridval process in this state can be defined as 

ki = max{Xti^j \ 5s.f, j = 1, 2, ...i - 1 }, (3) 

in which ti = i * 5s.f (i € {1, 2, ...t]}), and Atij = ti — tj. 
After Thar + Wrar time, a failed UE in any state will chose 
a number between [1 — Wbo] with equal probability So 
we can represent the probability of transmission from state j 

if a > 1 & Bo_maXj > ti 
if a > 1 & ti_i < Bo_ maxj < ti 
if a < 1 & Bo_ maxj > ti 
otherwise 

(4) 

In which the parameters related to each state is defined as 

a = Atij \ 5s.f, j = 1 , 2, ...i - 1 , 

Bo_mi7ij = tj + Trar + Wrar, 

Bo_maXj = Bo_mi7ij + Wbo, 

AWbo = {Bo_maxj — Bo_mi7ij) + 1 . 

We will explain the usage of transmission probability in an 
example in this section. 

As it comes from Eq. loop back or backward steps in 
this K-Mc model is not possible. Having Pij function, we 
can define the transition probability matrix of (p * 77 ) size as 


to the state i as following 


^s.f 

AWbo ’ 
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PiA 

Pi,2 

^4.3 

0 . 
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■ (5) 
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Where 2 < i < 7 } and 1 < J < * - 1. Having the Pij matrix 
we can estimate the expected number of arrivals at time t in 
state as defined in Eq. O where, 

ki 

z['^\t, h) = E * zf-^\t, %), 2<7i<W, (6) 

j=i-l 

W ki 

^ Zi(t, ki) = zl^\t) + E E * zi;-^\t, %). (7) 

n—2j—i—1 

In which Z^'^\t,kj) is the expected number of failed UEs 
arrival from state j who are attempting for re-transmission 
in current state i, and their population is defined as 

zj;) [t, k,) = zf [t, k, )*Pj„ 2 < 71 <w. (8) 

Pjj is the failing probability in state j. Let us define success 
probability Pi^ as the probability that exactly each preamble 
mi in ith{l < i < tj) state is chosen by one user, given the 
expected number of arrivals, Zi, as 

P^. = P{mi = l\Xiit, k) = Zi, mi e {1, 2...M}}, (9) 

so we can write 

Pij, = P{mi > l\Xi{t, k) = Zi, mi € {1, 2...M}} = l-Pi,. 

( 10 ) 

Perhaps expanding Eq. O worth to practice before we start 
deriving the other parameters. As it comes from Eq. (| 6 ]l the 
T^th ^jval at time t in current state is depend on failed 
population UEs from (n — 1)*^ time preamble attempt in 
previous states. Eor example, if the system at time t is in 
state z = 4, given 5s.f = 10ms, AIUbo = 21ms, from 
Eq. 0 we can compute the number of associated past states 
^4 = 3 and their transition probability P 4 J, {j = 1, ...3) as: 
F4.1 = 0.1, P 4^2 = 0.5, P4_3 = 0.4. Using Eq. (| 6 ll, O the 
expected number of arrivals in this state which is in associated 
with the new arrival (t) and failed users in previous states 
1, 2 & 3 can be denoted as 

Z4(f, ki) = zf V) + zi^\t, ki) + zi^\t, ki) + zf\t, ki), 

( 11 ) 

in which, 

zf\t) = Zi, (12) 

Zf\t, ki) = P4.1 * Z^'^it, ki) + Pi,2 * zi^^it, k2) 

+ Pi,3 * Z^^(t, fca), 

z!t\t, ki) = Pi,2 * Z^f(t, k2) + P4.3 * Zf^{t, ks), 

Zi^\t,ki)=Pi, 3 *Zifit,k 3 ). 

Where, Z4 is the contribution of initially distributed users in 
4*^ state at time t. Note that the rest of higher order failed 


















TABLE I. 


Parameters 


where. 


Parameter 

Value 

M (total number of preambles/RA-TS) 

54 

Ti?Ai?(Msg 2 response delay) 

2 ms 

^RAR (Msg 2 response window) 

5 ms 

6 s. f (RACH sub frame interval) 

10 ms 

if PRACH Configuration Index = 6 

W (maximum transmissions) 

10 

^BO (Backoff window) 

20 ms 

T (latest duration of the users activation) 

lOsec 

77 (last Markov state / RA-TS ) 

T 

^s.f 

H 2 H and M 2 M call probability distribution function 

Typel & Type 2 [IJ 


UEs (n > 4) are not possible in this state. As it comes from 
the above derivation we can easily compute the other states 
population recursively. 


B. Performance Metrics Analysis in RACH Procedure 

In this subsection we analyze the performance metrics in 
RACH process which is in associated with two types of users 
H2H and M2M. To begin with, we start with clear definition of 
performance metrics specified in 3GPP TR 37.868 as follows 

m. 


h) = ^ h) + zf h). 

n—1 

(18) 

Collision probability, defined as the ratio between 
the number of collided users and the total number 
of RAOs (with or without access attempts). If two 
(or more) users select the same preamble at the same 
RA-TS, the eNB will not be able to decode any of 
the preamble requests; hence, eNB will not send RAR 
msg. 

Having RAOs = rj * M, by using Eq. (fTOl l and Eq. 
dnii-® we can derive the collision probability as 


Pf = 


Eti EZi it, h) + (f, h) 


r] * M 

in which 

z. ^ exp{ - — - 


(19) 

( 20 ) 


)), 


1) Access success probability, defined as the ratio be¬ 
tween the number of successful users who complete 
the RA procedure within W number of preamble 
transmissions and the number of active users. Which 
should not be mistaken by the probability of success 
in respect to the RAOs (random access opportunity) 
as worked in a, M- It can be written as 

„ No. of success UEs within W retry 

Pg = -. (13) 

Total no. of active UEs 

In this work we have assumed the success preamble 
in step2 as a successful user. Eor the corresponding 
number of contending MTC, and hu¬ 

man type call arrivals, ZP'^^{t, ki), at time t in state 
i, the number of success UEs can be defined as 


Z^^it, h) = ""(f, h). 


(14) 

In which, the probability of success for H2H and 
M2M communication in state i is defined as a 


H 2 H 


{rrii = ki) = E[Zi^ (f, ki)] 








(15) 
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M 2 M 




m, e { 1 , 2 ...M}} 

_ /-r_ 

^ m’ 




(16) 


Therefore, the access success probability Eq. (HI can 
be computed as 


E r) 7 ^ 

Z^n=l 


Pg = 


rH2H{n 


\t,ki) 


ZM2M(n)^^^ki) 


EUZ,{t,ki) 


(17) 


3) 


4) 


( 21 ) 

Statistics of number of preamble transmissions, 

defined as the CDE of the number of successful users 
who connect to the eNB within r {r < W) number 
of preamble retransmission to the total number of 
successful users. Which can be written as 


FpiP <r) = 




( 22 ) 


Access success delay, defined as the expected delay 
for successful users who complete the RA procedure 
within W number of retransmission preamble to 
the total number of successful users. Which can be 
denoted as 


Z^i=l /-^n=l 

E ri y{' 

i=l Z^rt=l 


E[t] = 


rin) 


it,ki) 


(23) 


in which, ttI"! is the average success delay at 
nth retjansmission in state i. Basically users delay 
comprises of three terms 
User’s delay = delay 

proc. 

delay6„^„„g -f delay 

arture‘ 

Without considering the factor of barring delay, let’s 
consider those UEs who failed in previous ki states 
and their preamble transmission is successful in 
state i. We can write the expected access success 
delay for successful users in state as 


‘I - >1 ^ 5 


(24) 


where the expected success access delay for new 
arrivals ( 77 *-^^ delay of departure), in state i can be 
estimated as expected number of success UEs in this 
state multiplied by their average delay for departing 
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Fig. 3. Scheme of JA & DA splitting preamble methods in RACH procedure. 


between i — 1 and state, (Ss.f/2), and delay of 
receiving RAR message, (Trar + Wrar) as 

* +Trar + Wrar), (25) 

and the average access delay for higher order pream¬ 
ble retransmissions in this state can be estimated as 

j=i-l 

* {kj * 5s.f + Trar + Wrar)- (26) 

In which, term kj*5s.f is an extra delay for those UEs 
who collided in previous state (ki < j < i — 1) 
and re-attempt in state i. 

5) Statistics of access delay, dehned as the CDF of 
the delay for successful users who complete the RA 
procedure within w number of retransmission time 
divided by total number of successful users. Which 
can be denoted as 


FV(t < w) = 


L^i—1 2-^n—l ^ 


in) 




(27) 


IV. Splitting RACH Resources For M2M And H2H 
Communication And Parameter Optimization For 
Proposed JA Scheme 

According to H] when MTC and H2H calls share the 
RACH resources, they experience the same access collision 
probability. Providing separate RACH resources for the H2H 
and MTC devices (DA method) can reduce the number of 
collisions, however as we discussed in IS JA method has 
better performance in reducing the collision rate for MTC calls 
in comparing to the DA algorithm. In this section we briefly 
describe the mentioned preamble splitting RAN-overload con¬ 
trol algorithms, and then we develop the access success delay 
performance metrics through analysis. The target in this section 
is to And an optimum number of preamble allocation to H2H 
calls, x\, (1 < z < p), in order to maintain the required 
access success delay for machine type applications below a 
given threshold using JA method. 


A. DA and JA RAN-Overload Control Methods 

As it’s depicted in Figure [3 each state z contains M 
preambles M = {toi, m 2 , ...M} shared between human and 
machine type devices. In DA algorithm, users have access 
to a fixed number of allocated preambles and contend. For 
example, in this figure x = a and M — a is the number 
of preambles allocated to the H2H and M2M communication 
respectively. Whereas, in JA algorithm, H2H users have access 
to X number of preambles [x is a random variable) and M2M 
users contend with them using the access permission to the 
whole M preamble set. Which means that, the x number of 
preambles will be shared between M2M and H2H, whereas 
still these are M2M users who have access to the M — x of 
them. Let Ai be the average number of H2H new arrivals in 
state z, g{t) be the probability distribution function of RACH 
requests generated by the M 2 M calls, and Nmtc be the total 
number of M2M active users during [0 — T] time. The number 
of M2M new arrivals sending their preamble request in state 
z can be defined as defined in m. 

Ni = Nmtc g{i) dt, (28) 

We can write the expected number of new arrivals in z*^ 
RA-TS (z S {1,2, ...z;}) with respect to DA and JA method 
respectively as follows 

Z_DA[^\t)=X,\a + N,\M-a, (29) 

in which Ai and Ni users contend for a and M — a number 
of preambles individually, and 

Z_JA\^\t) = (Ai -I- — )\m-x^ (30) 

where (Ai + Ni^) number of initial users share x number of 
preambles together. Therefore, we can re-write the population 
of z*^ state in general (either for DA or JA method) as follows. 

Z_DA[JA]i{t, ki) = Z_DA\JA]f\t)+ (31) 

ki W 

j—i—1 n—2 

One should notice that we can And the expected number of 
arrivals by considering this fact that even though the new ar¬ 
rival of M2M users follow a Beta distribution, yet their higher 
orders in back-off procedure follows the uniform distribution. 
Therefore, this is easy to calculate the performance metrics 
with respect to the DA/JA method, that we neglect to write 
here. 


B. Preamble Allocation Optimization for M2M Calls Given 
the Access Success Delay Threshold & Number of Preambles 
per RA-TS 

An optimum value for x], that minimizes the access 
success delay for M2M calls using JA algorithm for a given 
maximum threshold delay that can be denoted as 

min{E[T^^^ {xj ,M)]} with respect to Xi, 


(32) 
















where, 
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in which from Eq. (l26l l. we can write 

^[p] = ( o -^) * i^s./+ + VEfiAfi (34) 


and, 
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Fig. 4. Estimated number of M2M UEs arrival after retransmission 
( 2 -M 2 M(n)) aiTival simualtion (fF = 10). 


Simulation Vs Analysis, lambda=10, T_PER=10ms 


Eq. ( |3^ is too complex to use. Let’s denote optimum xl such 
that 


where the two convext functions A, B are the numerator and 
denominator in Eq. (l3^ . Using the Jensen’s inequality for any 
$ convex function and random variable X we have 

$(E[X]) < E[$(X)], (36) 

let’s assum the expected number of M2M and H2H new 
arrivals in state i is Ni and Xi respectively. So we can write 

A,ix,,M)> (37) 

a,(x.)(l - + 7,(x,)(l - 

Xi M — Xi 

and 
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Fig. 5. Estimated number of H2H UEs anival after retransmission 
(^/f 2 /f{n)) aiTival vs. simualtion (W = 10). 


argmin{E[Tjf^^(4, M)]} = rmn{E[Ti]A{xl M), B-\xl M)} 


Bi{xi,M) > Ai{xi,M) + Xiexp{^^). (38) 

Xi 


Where ai{xi) = Xi + Nij^, ji^Xi) = TV, 


M—Xi 


. The R.H.S 


j . ’*M’ M 

of the Ea.lBTb is the minimum expected value for function A. 
Without loosing the generality, by finding the minimum value 
for the function A and maximum value for function B we 
can minimize the whole fraction. To maximize the function B 
we can let the system to reach its steady state in which the 
number of arrivals per RA-TS will reach to W times of its 
new arrivals. Therefore, by numerically solving the following 
equation we can obtain xj with respect to (j) as 

solvet ’>X.SlA(.hM,K.N.)l ^ , 

R.H.S[B{xl;M,W *Xi,W *N,) 


V. Numerical Results 

In this section, we aim to show the accuracy of RACH 
performance metrics approximation using the proposed K-Mc 
model in comparison to the simulation results. The analysis 
results include the optimal number of preamble allocation 
to H2H calls (i.e., x^) in PRACH, in order to avoid excess 
success access delay for machine type applications. Under 
the proposed JA method, we compare the analytical results 
with simulation to see for how many Nmtc of simultane¬ 
ous machine devices in a cell, the system can support the 
access delay requirements for M2M calls. Eigures dUl, © 
demonstrate the number of M2M and H2H arrival vs. their 
new arrival over RA-TSs. Where we compared the analysis 
results (i.e., K-Mc model Eq.©) and simulation. As it comes 
from these figures, the actual number of active M2M/H2H 




























Fig. 6. M2M access success probability with DA method with 

><h2h = 0.5. 
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users is very well matched with the analytical model. From 
Fig. @ one can notice that although M2M new arrivals is 
modeled with a Beta distribution, however, since the back¬ 
off procedure is a uniform distribution over RACH slots, the 
number of M2M arrivals after W times retransmission is a 
linear product of W in steady state. This phenomena 
is also applicable on H2H users with Poisson arrival, where 
A = 10 is the average rate of FI2H new arrivals, which reaches 
to A * FF = 100 per RA-TS in steady state. 

In Fig. (O and © we represent the results of access success 
probability for M2M users, with respect to DA and 

JA method respectively, where the average rate of FI2H calls 
per second is A = 0.5. We have verified the same trend for 
pM 2 M Q j < A < 1. As it comes from the Fig. (l6|, 

the probability of access success degrades with respect to the 
number of preamble allocation to H2H users, x. As x increases, 
M — X lesser number of preambles will be allocated to the 
M2M users. Where in JA method (Fig. ©), the success access 
probability falls down with a better slope. The reason behind 
this result is that in JA algorithm the number of x preambles 
is being shared with H2H and M2M users to content. As the 
result the M2M users can get success in lesser number of 
preamble transmission with JA algorithm than in DA (please 
see Fig. (l8]l). Following the results from the access success 
probabilities, we can expect a better access success delay in 
JA method in comparison to the DA method. Which, as it 
is depicted in Fig. the results of access success delay 
with respect to JA method is a smooth convex function with 
regards to x. Where, in Fig. (fTOl i the success delay for M2M 
users is dramatically increasing by allocating more number 
of preambles to the H2H users. For computing the success 
access delay for users, one should note that this value is 
composed of the delay of departure and delay in receiving 
success RAR acknowledgement for Msg2 from eNB in RA 
procedure. In this paper, we have assumed that the delay of 
other handshaking messages (Msg.3-4) in RACH for allocating 
the physical uplink channels is constant for all type of users. 
Therefore, the required success access delay threshold, </>, for 
M2M users is considered for successfully receiving the RAR 
message after sending the preamble request (Msg.2). 

The results of optimum preamble allocation to H2H users. 


Fig. 7. M2M access success probability with JA method with 

^h2h = 0.5. 



Fig. 8. CDF of success preamble trasnmission; Xh 2 h = 0.5. 


x'^, with respect to ^ is shown in Fig. As we can 

deduce from this figure, the optimum number of allocating 
preambles to H2H users per RA-TS decreases with lesser 
access delay threshold requirement, (j). For example, if there are 
Nmtc = 5000 number of machine devices in a cell have delay 
requirement of </> < 20ms, only x = 5 number of preambles 
out of M available preambles in a RA-TS should be allocated 
to H2H and be shared with M2M subscribers. If we consider 
a real scenario in which there are Nmtc = 30000 machine 
type devices with 100ms success access delay, which requires 
receiving the RAR Msg. in less than 40mfl there will be only 
X = 3 preambles to allocate to the H2H users. Under the bulky 
arrival of M2M users Nmtc = 30000 with average arrival rate 
of H2H users 0.1 < A < 1, M2M users can get success access 
to the eNB only after 35ms. 


^Considering the delay of receiving Msg.2 consumes 1/3 of total access 
success delay in RA procedure. 
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Fig. 9. Access success delay for M2M calls using JA method, 


Fig. 11. Optimum no. of preamble allocation to FI2Ft users, (a;t ((/>)) 

, with respect to the given M2M access success threshold (p. 



Fig. 10. Access success delay for M2M calls using DA method, 


VI. Conclusion 

In this work, we have proposed a closed-form Marko¬ 
vian model for splitting preamble allocation RAN-overload 
algorithms JA and DA methods. Under the same model we 
have derived a closed-form analytical expression for computing 
performance metrics in random access procedure including 
collision probability, access success probability, access success 
delay, CDF of preamble transmission and CDF of access 
success delay for H2H and M2M users. In addition, we 
examined the JA method for reducing the access success delay 
for M2M users in the random access procedure. Under the JA 
algorithm we have suggested a formula for the optimal number 
of preamble allocation to the H2H subscribers in PRACH. 
Finally, we have shown that the numerical results from the 
proposed order Markovian model verifies the simulation 


results. 
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